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Traditionally quantum tunneling in a static SQUID is studied on the basis of a classical trajectory 
in imaginary time under a two-dimensional potential barrier. The trajectory connects a potential 
well and an outer region crossing their borders in perpendicular directions. In contrast to that 
main-path mechanism, a wide set of trajectories with components tangent to the border of the 
well can constitute an alternative mechanism of multi-path tunneling. The phenomenon is essen- 
tially non-one-dimensional. Continuously distributed paths under the barrier result in enhancement 
of tunneling probability. A type of tunneling mechanism (main-path or multi-path) depends on 
character of a state in the potential well prior to tunneling. 
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I. INTRODUCTION 

Despite phases in Josephson junctions [l| are treated 
as macroscopic degrees of freedom they can exhibit quan- 
tum properties 0. In particular, quantum tunneling of 
those variables across a potential barrier is possible ■ 
Tunneling in a single Josephson junction is similar to a 
conventional one-dimensional quantum mechanical pro- 
cess. In this case the tunneling mechanism is described 
by theory of Wentzel, Kramers, and Brillouin (WKB) 7]. 
Tunneling occurs from a classically allowed region vifhich 
is a conventional potential well where energy levels are 
quantized 8-11]. Quantum coherence between potential 
wells was demonstrated [l2| - [T5| . 

Besides single Josephson junctions superconducting 
quantum interference devices (SQUID) are also a mat- 
ter of active investigation for many years [l6l - l23| . A 
SQUID consists of two Josephson junctions and, there- 
fore, represents a two-dimensional system where macro- 
scopic quantum tunneling is also possible. Tunneling in 
multi-dimensional systems is well studied [23 - [28j . There 
is the certain underbarrier path (or a few paths) where a 
wave function is localized and it decays along the path. 
This main-path tunneling is described by a classical tra- 
jectory in imaginary time and it is generic with a con- 
ventional WKB mechanism. 

In contrast to that, a different scenario of tunneling 
in a static SQUID is possible. Instead of localization on 
a main path an underbarrier state is distributed over a 
continuous set of paths. Multi-path tunneling cannot be 
described in terms of a classical trajectory in imaginary 
time. 

A realization of a particular tunneling mechanism 
(main-path or multi-path) depends on type of a state 
in the potential well prior to tunneling. Suppose, in the 
classical language, a state in a well to have a momentum 
component orthogonal at some point to a border of the 
well (a normal reflection). The main path starts at that 
point and continues under the barrier. When a state in 
the well has also a component which is tangent to the bor- 
der at some point (a particle hits the border from some 



angle) a scenario is different. Around that point there a 
set of continuously distributed classical paths which go 
under the barrier. 

A goal of the paper is to propose the different as- 
pect of tunneling in a static SQUID which is multi- 
path mechanism. The phenomenon is essentially non- 
one-dimensional. As shown in the paper, a multi-path 
mechanism can result in a larger probability of tunneling 
compared to main-path. Different paths interfere and 
a method of classical trajectory in imaginary time, in 
contrast to main-path, is not valid. In experiments it 
is possible to determine a contribution of multi-path ef- 
fects to a total probability of tunneling. A symmetric dc 
SQUID in zero magnetic field and without dissipation is 
considered below. 

In Sec. In] a formulation of the problem is given. In 
Sec. mil general arguments are used to explain the phe- 
nomenon and to estimate the effect. It is shown that 
multi-path tunneling through a two-dimensional static 
barrier reminds photon-assisted tunneling across a non- 
stationary one-dimensional barrier. In Sec. IIVI an exact 
solution of the semiclassical problem is done with the 
use of a certain model coupling between the junctions in 
a SQUID. The exact calculations confirm the estimate of 
the effect performed on the basis of general arguments in 
Secimi 



II. FORMULATION OF THE PROBLEM 

We consider a dc SQUID, consisting of two identical 
Josephson junctions with phases ipi and ip2^ with no dis- 
sipation when the two junctions are inductively coupled. 
A classical behavior of phases corresponds to conserva- 
tion of the total energy 



En - 



Ej_ 



dt 



d(p2 
dt 



+ Ej 



cos ipi 



■C0S(/?2 - JfiVl + f2) + -^ifi - V2f 



(1) 



2 



Vo(y) 




Vo(x) 



Vo(y) 



FIG. 1: Uncoupled junctions, a = 0. The total particle en- 
ergy E is fixed. The maximal tunneling probability in the x 
direction relates to a minimal excitation of the y motion to 
provide a maximal energy in the x direction. 



where the Josephson energy Ej = hlc/2e, the plasma 
frequency uj = y^2eIc/hC, and the coupling parameter 
/3 = 2ttLIc/^o are introduced. Here Ic, L, and C are crit- 
ical current, inductance, and capacitance of each individ- 
ual junction. The magnetic flux quantum is $o — T^hc/e. 

Below we consider large f3 and the total current / close 
to its critical value 



1 



2Ir 



< 1 



and introduce new variables by the relations 

VI ^ ^ + {3x ~ 1)^2(1 - I/2Q + j 

V2 = ^ + {3y- l)V2(l-//2/e) + ^. 
Below time is measured in the unit of to defined by 
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The energy ([T]) takes the form 
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The potential energy is 
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where Vq (x) = x^ — x^ and the coupling parameter is 
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B in Eq. (|6|) is called semiclassical parameter. When B 
is large the phase dynamics is mainly classical. Below we 
consider that case, 1 <C i?. 




FIG. 2: Uncoupled junctions, a = 0. A part of the total 
energy E is given up to an excitation of the y motion resulting 
in a reduced tunneling probability in the x direction compared 

to Fig. m 



III. TUNNELING IN TWO DIMENSIONS 

A classical dynamics of phases in a SQUID is described 
by Eqs. ([5]) and ([7|). The effective particle moves in 
the classically allowed region, in a vicinity of the point 
X = y = 0, which is restricted by the potential barrier. As 
known, the particle can tunnel through the barrier result- 
ing in experimentally observable phase jumps. Character 
of tunneling depends on coupling strength a between the 
two junctions. When the coupling is strong, 1 ^ a, one 
can easily show that the last term in the energy ([1]) dom- 
inates Therefore in this case the particle tunnels 
along the direction x = y since fluctuations around this 
path cost a large energy. In contrast, at a small cou- 
pling, a ^ 1, the junctions are almost independent and 
there are two different tunneling paths, along the x or y 
direction [l7l |. 

Therefore at a very strong and a very weak coupling 
between two junctions the system behaves as effectively 
one-dimensional. Essential features of two dimensions are 
exhibited for an intermediate coupling, a 1. Below we 
study a formation of that regime starting from a region 
of small a. 



A. Uncoupled junctions 

Let us consider first zero coupling between two junc- 
tions (a = 0). In this case Eqs. dSJ and (O describe two 
independent particles in one-dimension potentials shown 
in Fig. [TJ When the semiclassical parameter B is large 
the potential barriers in Fig. [T] are hardly transparent 
and a number of discrete levels in the wells is large. 

Suppose tunneling to occur in the x direction. We 
introduce the dimensionless energy E 
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The total particle energy E in Fig. [T] is a sum of ones cor- 
responding to motions in the x and y directions. A maxi- 
mal tunneling probability of a particle, with a fixed total 
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FIG. 3: At 5E = the main path (arrowed curve) con- 
nects two classically allowed domains, the potential well and 
the outer region. The thick curves relate to the condition 
V{x,y) = E. A deviation from the main path (dashed curve) 
does not lead to the well. On the line c the y component of a 
force is zero, dV/dy — 0. 



energy E, is realized when the energy of the x motion 
has a maximal possible value. This situation is shown in 
Fig. [T] In this case a motion in the direction perpendicu- 
lar to tunneling (the y direction) is not excited pertaining 
to a lowest level in y. In the classical language, the par- 
ticle hits a border of the well, V{x,y) — E, with zero 
tangent velocity. 

For comparison, in Fig. [5] the total energy E is dis- 
tributed in a way that the y motion is excited. In the 
classical language, the particle hits a border of the well 
with a finite tangent velocity. This results in less prob- 
able tunneling in the x direction since tunneling occurs 
from a lower level, E — SE, as shown in Fig. [2] The part 
SE of the total energy relates to the tangent motion. 



B. Weakly coupled junctions 

Suppose the junctions are weakly coupled so that 
a <^ 1. In this case a motion in the total potential ([7]) 
cannot be reduced to two independent ones as in Figs. [T] 
andlH Now the entire system of levels in the total poten- 
tial (I7|) accumulates the levels of the x and y channels in 
Fig. [TJ If there are five discrete levels in each well, Vo{x) 
and Vb(y), then the potential V{x,y) contains 25 levels. 

Let us analyze specificity of two-dimensionality in tun- 
neling. First, the coupling constant a is partly accounted 
for in the parameter B ^ {1 + a)^/^ in Eq. ([5]). At small 
a this leads to a linear in a reduction of tunneling prob- 
ability. The last term in the potential ([T]) is proportional 
to since y ^ a. So we calculate an correction to 
the tunneling probability. It is small at a <C 1 but it 
becomes significant at a ~ 1. 

Below in this subsection we propose some non-rigorous 
arguments which help to understand what happens under 
the barrier in two dimensions. 

When a particle hits a border of the well with zero tan- 
gent momentum, 5E ~ 0, an underbarrier wave function 
is localized on the certain classical trajector y w hich is or- 
thogonal to the borders of classical regions [24| - |28| . This 
trajectory is shown in Fig. |3]by the arrowed curve. The 
trajectory is driven by the y force directed away form the 



FIG. 4: Underbarrier paths at 5E 7^ connecting two classi- 
cally allowed domains as in Fig.O A finite tangent component 
results in the multi-path region close to the well which goes 
over into the subsequent path indicated by the arrowed curve. 
The main path of Fig. [3] is shown by the dashed curve. 



line y — ax /{I + a) where it changes sign. That line is 
marked in Fig. [3] as c. Not far from the outer region the 
X force, dV/dx, attracts the particle to the line c as it is 
tilted. 

A wave function is proportional to exp(iS'/fi,) where S 
is a classical action [7]. According to Maupertuis' prin- 
ciple dH, along the classical trajectory one should put 
S = Sci where 



Sci = ihB / dl^V{x,y)-E 



(10) 



and dl is an element of the classical trajectory. This 
mechanism can be called main path tuneling since an- 
other path does not lead to the well in Fig. [3] deviating 
from it as the classical trajectory shown by the dashed 
curve 0, 113] ■ Eq. ([T0| is generic with a conventional 
WKB result if to consider it along a curve but not a 
straight line. 

A scenario changes when a particle hits the border 
of the well with non-zero tangent momentum, 6E ^ 0. 
Then an individual trajectory, related to the regime pH)) 
and shown by the arrowed curve in Fig. 21 is impossible at 
least in a vicinity of the well where a tangent component 
of the momentum is finite 
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dy 

and the imaginary momentum in the x direction is 
diS 
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dx 
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In this case an attempt to adjust a classical trajec- 
tory to a tangent component at the well border results 
in a deviating curve similar to the dashed one in Fig. [3l 
Therefore an underbarrier density cannot be localized on 
a particular classical trajectory along the entire under- 
barrier region. Close to the well that trajectory goes 
over into a wide set of paths as shown schematically in 
Fig. 131 According to Feynman [sO], the state related to 
([TT|) and (fT2)) is a superposition of a wide set of classical 
paths. 

As follows from Eq. ((T^), a finite 5E reduces the wave 
function in a vicinity of the well. In contrast, far from 
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the well the wave function exp{iSci/h) on the arrowed 
curve in Fig. U] is larger compared to the dashed curve 
{6E = 0) since the former is closer to the line c where 
V{x, y) is smaller. In other words, the arrowed curve in 
Fig.|4]is shifted to a more transparent part of the barrier. 

When the energy E is close to the bottom of the well 
the fraction of the trajectory above the curve c in Fig. [3] 
is small because the well shrinks. 

As follows from above, there are two opposite effects 
on tunneling when the total energy in the well is redis- 
tributed between a normal {E — 5E) and a tangent {SE) 
motions. 

(i) The reduction of tunneling is connected to a vicin- 
ity of well where the wave function is suppressed due to 
sinking of an energy level down to E — SE. As follows 
from Eq. , the exponential reduction of the particle 
density is exp{—ciB5E) where ci is a positive parameter. 

(ii) The enhancement of tunneling is due to a shift of 
the subsequent trajectory to a more transparent part of 
the barrier. A reduction of ^(a;, y) in that region is pro- 
portional to axy. Since a typical a; ~ 1 and y ~ V SE, 
the exponential enhancement of the particle density, ac- 
cording to Eq. ([To]), can be estimated as exp(2c2Ba^/6E) 
where C2 is a positive parameter. 

A total effect on the tunneling probability T{SE) is 
defined by a product of the two exponents 



T{SE) - r(0)exp(-ciB(5£;) exp{2c2Ba\ 



(13) 



At SE ^ o? the expression (jl3p reaches a maximum, 
r(0) exp(i?a^C2/ci), which manifests an exponential en- 
hancement of tunneling due to a finite SE. The param- 
eters Ci and C2 are approximately of the order of unity. 
Whereas a is small the semiclassical combination Bo? is 
large. 

It is amazing, that the above conclusions, drawn on 
the basis of general arguments, are confirmed (excepting 
some details) by exact calculation in Sec. IIVI where the 
values of c\ and C2 are specified. 



C. Analogy with photon-assisted tunneling 

Multi-path tunneling through a two-dimensional static 
barrier reminds photon-assisted tunneling across a non- 
stationary one-dimensional barrier 31]. The latter also 
consists of two parts. The first one is an absorption of 
quanta with an exponentially small probability analogous 
to the first exponent in Eq. ((T5)) . The second one is tun- 
neling in a more transparent part of the barrier (a higher 
energy) with an enhanced probability corresponding to 
the second exponent in Eq. p^ . A total probability is 
also determined by a maximum of the product. 

In the both cases a particle finds a more transparent 
part of a barrier being initially pushed either by a tangent 
motion (at the end of the multi-path region in Fig. 2]) or 
by quanta absorption. 



IV. UNDERBARRIER WAVE FUNCTION 

To quantitatively study the problem of two- 
dimensional tunneling one should solve the Schrodinger 
equation with the exact potential ([7]). Since the potential 
barrier is almost classical one can apply a semiclassical 
method. With an exponential accuracy the wave function 
is 



V' ~ exp(ii?o-\/2), 



(14) 



where the exponent is a large classical action measured 
in the units of Planck's constant. As follows from the 
form ([5]), a{x,y) satisfies the Hamilton- Jacobi equation 

in 
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where f{x) — 2ax/{l + a). When a = the variables in 
Eq. ([T5|) are separated and a solution can be easily ob- 
tained. In our case there a substantial cross-term f{x)y 
in Eq. (1151) which mixes the modes. 

We consider a small coupling a when a deviation of 
the variable y from the tunneling path is small. For con- 
venience, below a model coupling between the variables 
X and y is introduced. Namely, instead of the linear in x 
function we use 
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where xq is chosen between 02 and ai in Fig. [T] The 
choice of the form (fT6l) does not contradict to main ar- 
guments of Sec. 



A. Hamilton-Jacobi approach 

At a small coupling, a ^ 1, we consider a transition 
through the barrier along the direction x. The transition 
occurs with a small deviation of y from zero position. 
In this case one can put Vo{y) = in Eq. (fT5|) . The 
classically allowed motion in the well Vb(x) + y"^ < E is 
described by a real action a{x,y). The particle probes the 
potential in the classically allowed region where motions 
in X and y directions are independent. Therefore the level 
quantization is determined by two one-dimensional wells 
and the action is a sum of two one-dimensional parts 
defined by a solution of Eq. ((T5| with the condition (fT6)) 
at X < Xo . A continuation of this action from the well to 
under the barrier at x < xo results in 



a{x,y)^i dxiy^Voixi)-E + SE+ / dyiJSE-yl 

Jxa JO 

(17) 

The part SE of the total energy E relates to the tangent 
(perpendicular to the tuneling direction) motion. The 
fraction SE is determined by a state in the well from 
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FIG. 5: Real momenta are indicated by arrows. A maximum 
of the particle density under the barrier is reached on the thick 
curves y{x). (a) Tunneling from a state with zero tangent 
momentum, SE — 0. (b) Tunneling from a state with non- 
zero tangent momentum, 5E 7^ when yo = V SE is finite. 
Real momenta at the well {x < 02) are shown on a classical 
trajectory only. 



which tunneling occurs. At a large B the energies E and 
SE are almost continuous with a discreetness 1/B. 

A solution at xq < X can be found by the method of 
variation of constants (291. It has the form 



a{x,y)^i dxiy/Vo{xi) - E - a'^ + e{x,y) (18) 

Jxo 

rv 

+ / dyi^ e{x, y) - {yi - of + iF{e{x, y)\, xn < x. 



For a given function F{e) the function e{x,y) is deter- 
mined by the condition da/de = 
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(19) 



a ~ {yi ~ ay 



which is independence of a on "constant" e{x,y). Now 
derivatives of the action have the simple forms 

= t^Voix) -E-a^+ e{x, y) (20) 

da{x,y) 



dy 



= \/e(a;,y) - [y - af 



According to matching of da{xQ, y)/dy given by (fT7|) and 
([T5|. e(xn.v) ~ 6E + a^ — 2ay. This can be used to deter- 
mined the function F(e) if to express y through e{xo,y) 
and to insert that into Eq. (fT9|). As a result, the function 
e{x, y) is defined by the condition 



dxi 
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B. Probability of tunneling 

Not very close to the classical exit point a\ , determined 
by the condition V^(x) ~ E, one can write the left hand 
side of Eq. (I^TI) as iT[x) where 

dx\ 



t{x) 



(22) 



.0 ^JV^{x^)-E 

Now a solution of Eq. (PT|) can be easily obtained. At 

Xq <x 



y/{y — a)2 — e ^ y (y + 2Q!sinh^ T/2y — SE — a sinhr. 

(23) 

There is a remarkable underbarrier path, y{x), where 
da{x, y)/dy = 0. According to Eq. (PD|) . at xq < x 



y{x) = y sinh t{x) + SE - Q![coshT(a;) - 1]. (24) 

The path y{x) is shown in Fig. [5] by thick curves. 

The region x < a2 pertains to the classical motion in 
the well. At a2 < 2^ < a particle density decays away 
from the region —6E < y < 6E. At xq < x the under- 
barrier density reaches a maximum (with respect to y) 
on the curve y{x). At a finite 6E there are real momenta 
under the barrier indicated by arrows in Fig. [Sfb). At 
Xq < X at the region with real momenta the wave func- 
tion is exponentially smaller than on the curve y{x). So 
at Xq < X one can account for the path y{x) only. It is 
generic with the arrowed curve in Fig. 21 At oi < x the 
particle escapes from under the barrier. 

Below it is convenient to introduce 

fi{x)^e[x,y{x)]. (25) 
As follows from Eqs. ^ and ([24| . 

a/ fi{x) = a cosh t(x) — \I oP- sinh^ t(x) + 8E. (26) 

We define probability of tunneling with a fixed energy 
E and an exchange iE as 



' exp (^2iBV2 J 



^-(02,0) 

da[x,yix)] 



dx 



dx 



(27) 



The definition (f27| is analogous to Eq. (|13p . It is usefull 
to consider the probability of tunneling r(0) at 5E — 
0. To obtain it one has to expand the root in the first 
equation ([20]) with respect to e — and to substitute 
into Eq. With the use of Eqs. ^ and ^ we 

obtain 



r(0) = TwKB exp 



Ba^ 



2T(ai) - H- e 



-2r(ai) 



where 

^WKB 



exp 



-2BV2 / dxy^Va{x)-E 

J a-> 



(28) 
(29) 



is a one-dimensional WKB like expression. 
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FIG. 6: Modulus of the wave function at SE — o? . In this 
case y(x) = a at all a2 < a; < ai. 



C. Tunneling from a state with non-zero tangent 
momentum, &E 7^ 

We consider weak tangent momenta related to the con- 
dition i5i? ~ o? . 

According to Eqs. dU]), (HOI), and (gH), at 6E < 
the tunneling probability is 



V(bE) = r(0) exp ( -BbE^pi\T{a2)\ 



BSE , 
^'''SE 

(30) 

The first term in the exponent comes from a conventional 
WKB reduction of tunneling probability at x < xo when 
the energy is reduced by SE. The second term in the 
exponent is due to the region xq < x where the trajec- 
tory goes in a more transparent part of the barrier. The 
second term dominates at SE ^ and increases the 
tunneling probability. 

As follows from Eq. (1^51) . at SE close to a^, ii{x) = 



{a — vSE) / cosh t{x). This means that the increase 
of the tunneling probability, associated with the region 
Xq < X, is more effective at SE = o? . According to 
Sec. mil when E is not far from the well bottom the 
crossover between two regimes is shifted toward the point 
X = a-i. For this reason, below we choose the parameter 
Xq in Eq. to be close to a-i for simplicity. In this 
case, at SE close to a^, the tunneling probability is 



Y{SE) = r(0) exp 



1 



= -2r(ai) 



— 5- (a - V^)^tanhr(ai) 



(31) 



There is another expression of the probability which 
follows from Eqs. ([23) and 



V{SE) = EvKKS exp \Bo?^Pi 

[a-^fSEf , , ; 
5 tanh T\a\) 



T(ai) 



(32) 



wave function is plotted in Fig. [6] where y(x) = a for all 
underbarrier track, 02 < a; < ai. There are two branches 
under the barrier and an outgoing wave after the exit 
indicated by the dashed line in Fig. [HI 



D. Temperature dependence of tunneling 
probability 

At a fixed particle energy E a tunneling probability is 
determined by Eqs. ([32]) and (HH) with SE = a^. At 
a fixed temperature all energies contribute to tunnel- 
ing with Gibbs factors exp(— i?o/r). In the semiclassical 
limit, 1 <C S, there are many levels in the well. Therefore 
one can treat E as a, continuous variable and to optimize 
the probability r((5i?) exp(— i?o/r)) with respect to E 
accounting for Eq. This procedure gives the certain 
optimal energy Et from which tuneling occurs 
.27]. If to omit the term proportional to in V{SE) 
an optimal energy Et is determind by 



T(ai) 



ioTV2' 



(33) 



Imaginary time T{ai) is given by Eq. (j22[) where one 
should put now xq ~ a2. It follows that at T = the 
optimal energy is Et ~ 0. Et at the barrier top in 
Fig. [2] corresponds, at small a, to the critical tempera- 
ture To = h/i:to^/2 when the decay occurs solely due to 
thermal activation. 

Now one can express the tuneling probability F^ at a 
fixed temperature in the form 



Tt = F^°^ exp 



Ba' 
75 



1 — exp — 



T 



(34) 



(0) 



where the exponential coinsides with one of Eq. ([5T|) if 
to insert there SE — and the expression ([55]) . F, 
is a probability of a conventional tunneling when a tan- 
gent momentum in the well is zero. As follows from the 
expression (|34l) . F^' < Tt- This means that tunneling 
occurs from states with a finite tangent momentum. 

Since energy E has to be larger than SE{= a^) tem- 
perature cannot be too low, l/ln(l/Q!) < T/Tc < 1. 



E. Dependence of tunneling probability on current 

Let us consider first an almost continuous distribution 
of levels in the well (1 ^ B) as in subsection D. 

The coupling constant a is of the order of unity at 
I = Ir where 



2IcJ 



(35) 



The Gaussian form (15^ is analogous to Eq. (IT51) . The 
tunneling probability reaches a maximum at SE — as 
it has been predicted in Sec. IIIIl The modulus of the 



Since currents are close to the critical value 11^ the pa- 
rameter /3 is supposed to be large. A dependence of F7- 
on current is schematically shown in Fig. [T] The limit 
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FIG. 7: Dependence of tunneling probability at a fixed tem- 
perature Ft on current in the case of almost continuous dis- 
tribution of levels in the well (see insert). The dashed curve 
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FIG. 8: Dependence of tunneling probability at a fixed tem- 
perature Ft on current in the case of discrete levels in the 
well (see insert). The dashed curve is the same as in Fig. [T] 



corresponds to the conventional probability ln(l/r 



of a small a, considered above and defined by Eq. ([M]). 
corresponds to the left part of the curves in Fig. [T) Parts 
of the curves in Fig. [7] closer to 2/^ pertain to a large a. 
In this case tunneling occurs symmetrically, x —JU since 
a deviation from this path costs a lot of energy • 

The dashed curve in Fig. El is r^°^ It IS very easy 

to numerically calculate F^'' using a standard technique 
when a two-dimensional trajectory is normal to clas- 
sically allowed regions 0, H, [23l - [27j . F^^ depends on 
parameters (coupling strength, asymmetry of junctions, 
etc.) of a SQUID fabricated for measurements. For com- 
parison with experiments one has to know those values 

exactly. For this reason, we leave F^' as a schematic 
dashed curve in Fig. [7] planning its exact calculation in 
the nearest future for particular parameters of a SQUID 
with an intermediate effective coupling. 

The curves in Fig. [7] are taken at a fixed temperature 
which exceeds Tq = [2(1 - I / Ic)Y'^fiuj /2tt sufficiently 
close to 2Ic (large a) [17|. In other words, in a close 
vicinity of 2/c the decay is due to a thermal activation. 
We will consider details elsewhere. 

When the semiclassical parameter B is not too large 
a tunneling probability still remains exponentially small 
but a level distribution in the well becomes substantially 
discrete. This is a typical experimental situation [8-11]. 
The underbarrier process sets the certain energy 5E — c? 
which is given up to the tangent motion. According to 
Eq. (j31[) . the tunneling probability reaches a maximal 
value under that condition. The energy exchange is sim- 
ilar to one illustarted in Fig. [51 

On the other hand, the optimal bE may not be exactly 
fitted by discrete energy levels in the well which can be 
outside the Gaussian distribution (|3ip . This contrasts to 
a continuous distribution of levels at a large B when a 
proper level always exists. A mismatch between bE and 
the level structure results in reduction of the tunneling 
probability. 

For this reason, one can expect a "resonance" due to 
coincidence oi8E{= o?) with a distance between discrete 
levels. There are no sharp peaks but rather a shallow 
wavy structure in Vt illustrated in Fig. [51 Each minimum 



results from the "resonance" 
interval of currents 

bl 



and is smeared out at the 



(36) 



which follows from the Gaussian distribution in Eq. (|31|) . 
Under the condition p5l) a depth of each minimum in 
Fig. [51 is of the order of one. In Fig. [51 only two minima 
are shown. 

Strictly speaking, the semiclassical theory is not appli- 
cable to a case of moderate B and the above arguments 
are heuristic. 



V. DISCUSSIONS 

Quantum tunneling across a one-dimensional static po- 
tential barrier is described by WKB theory. Tunnel- 
ing through a multi-dimensional barrier is well studied 
[24 -[27 1 . Accordingly, the main contribution to a tunel- 
ing probability comes from the extreme path linking two 
classically allowed regions. The path is a classical trajec- 
tory with real coordinates which can be parametrized by 
imaginary time. The underbarrier trajectory is a solution 
of Newton's equation in imaginary time. The trajectory 
is given rise by a particle hitting normally (with zero 
tangent momentum) a border of the classically allowed 
region. In terms of discrete levels in the well tunneling 
occurs from a state with the analogous property. Under 
the barrier the probability density reaches a maximum at 
each point of the trajectory along the orthogonal direc- 
tion with respect to it. Therefore around the trajectory, 
which plays a role of a saddle point, quantum fluctua- 
tions are weak. The wave function, tracked along that 
trajectory under the barrier, exhibits an exponential de- 
cay generic with WKB behavior. The above mechanism, 
which can be called main-path tunneling, was explored 
for two-dimensional tunneling in a SQUID in Ref. jT^ . 

However, in some cases tunneling through multi- 
dimensional barriers occurs according to a different sce- 
nario which is far from being similar to WKB. An exam- 
ple of such situation in a SQUID is investigated in the 
paper. 
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When a state before tunnehng has a tangent momen- 
tum with respect to a border of the well there are no 
extreme points on it since the derivative of the wave 
function along the border is finite. This means that a 
tunneling probability is no more determined by the main 
underbarrier path but comes from a wide set of paths. 
Traditionally, a decay of a state with a tangent momen- 
tum is not considered since it does not correspond to a 
saddle point and, hence, the net contribution is expected 
to be averaged down to a small value due to mutual in- 
terference of trajectories. 

It is shown in the paper by general arguments and ex- 
act calculations that a non-zero tangent momentum does 
not result in reduction of tunneling probability. More- 
over, a related multi-path mechanism can exponentially 
enhance barrier penetration in a SQUID. In contrast to 
main-path mechanism, the multi-path tunneling cannot 
be calculated using a classical trajectory in imaginary 
time. Such a trajectory, connecting two physical (real) 
points, does not exist in that case. 

A possible way of interpretation of multi-path tunnel- 
ing is proposed in Sec. lIIIl In a vicinity of an enter under 
the barrier it is locally less transparent. A wide distribu- 
tion of paths in that region goes over into the single path 
which proceeds up to an end of the barrier. That path 
lies in a more transparent part of the barrier resulting in 
enhancement of tunneling probability. 

Multi-path tunneling through a two-dimensional static 
barrier reminds photon-assisted tunneling across a non- 
stationary one-dimensional barrier. In the both cases a 
particle finds a more transparent part of a barrier being 
initially pushed either by a tangent motion or by quanta 
absorption. 

The total energy E oi a state in the well is distributed 
as E — SE for a motion in the tunneling direction and as 
6E for tangent one. For a small coupling constant a the 
maximal tunneling probability corresponds to 5E ~ a^. 

An analysis of experimental data allows to define a 
contribution of multi-path effects by separation of the 
conventional probability F^-* from a total result. F^^ 
depends on parameters (coupling strength, asymmetry 
of junctions, etc.) of a particular SQUID fabricated for 
measurements. 

A SQUID should be with a large parameter /3. The cal- 



culations in the paper are done for a symmetric SQUID 
in zero magnetic field but this is not a principle restric- 
tion for observation of multi-path tunneling. A role of 
dissipation is to be studied. 

We see that in a multi-dimensional system (a SQUID 
is a two-dimensional example) tunneling mechanism can 
be different than one traditionally considered on the basis 
of a main trajectory. There are two necessary condition 
for unusual tunneling. First, tunneling should not start 
from a ground state, otherwise an energy exchange is 
impossible. Second, coordinates should not be almost 
separated when tunneling occurs practically along one of 
them (the limits of small and large coupling in a SQUID). 

One should expect a modified quantum tunneling at 
a finite temperature across a non-one-dimensional bar- 
rier when coordinates are not separated. In this case 
the conventional mechanism of a periodic trajectory in 
imaginary time with the period h/T is substituted by 
multi-path. 



VI. CONCLUSION 

Traditionally quantum tunneling in a static SQUID 
is studied on the basis of a classical trajectory in imagi- 
nary time under a two-dimensional potential barrier. The 
trajectory connects a potential well and an outer region 
crossing their borders in perpendicular directions. In 
contrast to that main-path mechanism, a wide set of tra- 
jectories with components tangent to the border of the 
well can constitute an alternative mechanism of multi- 
path tunneling. The phenomenon is essentially non-one- 
dimensional. Continuously distributed paths under the 
barrier result in enhancement of tunneling probability. A 
type of tunneling mechanism (main-path or multi-path) 
depends on character of a state in the potential well prior 
to tunneling. 
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